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G^ ■ Abstract 

CN I The different methods of reducing decoherence in quantum de- 

^ ■ vices are discussed from the unified point of view based on the en- 

. ergy conservation principle and the concept of forbidden transitions. 

[ Minimal decoherence model, "bang-bang" techniques, Zeno effect and 

^ ' decoherence-free subspaces and subsystems are studied as particular 

' examples. 

o : 
^ ■ 

1 Introduction 

Ch ■ 

■ The rapid development of quantum information theory, both in its theoretical 
^ • and experimental aspects [TJ |2] has posed new and renew old questions and 

challenges in the theory of quantum open systems jHlIll- Fulfilling extreme 
• ^ ! demands concerning the accuracy of controlled operations {unitary gates) 

^ I performed on microscopic quantum devices is crucial for any attempts to im- 

■ plement the ideas of quantum computation or generally quantum information 
processing at a reasonable scale. 

One can find in the literature different ideas of various technical complex- 
ity which should be useful for reducing errors due to decoherence processes in 
quantum systems interacting with an environment. We shall concentrate our 
attention on four of them : minimal decoherence model [3111; "bang-bang" 
techniques [7| , Zeno effect jHl IHl and decoherence-free subspaces and subsys- 
tems jlUt ITT] . The whole field of active quantum error correction methods 
|12j is not discussed here. An interesting question remains, whether these 
methods can be explained in a similar fashion also. 
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Mathematical relations between different methods of noise reduction in 
quantum systems have been already discussed in the literature ■ The main 
goal of this paper is to show that those different techniques can be under- 
stood from a unified physical point of view of two elementary principles - the 
energy conservation and the existence of forbidden transitions. Moreover, for 
the quantitative approach only the elementary time-dependent perturbation 
theory {Fermi Golden Rule) is used. One should mention, however, that for 
some cases more advanced and rigorous methods are known in the literature. 

1.1 Fermi Golden Rule 

Consider a quantum system with a Hamiltonian Ho perturbed by a "small" 
term V such that the total Hamiltonian reads 

H = Ho + V. (1) 

Assuming the spectral resolution of Hq in the form 

Ho\k >= Ek\k > (2) 

one can obtain the celebrated Fermi Golden Rule for the transition probability 
per unit time from the state |m > to the state \n > caused by the small 
perturbation: 

The formula (jHl), as stated, is rather singular and needs some additional expla- 
nation. In particular, taking sums over certain sets of initial and final states 
is usually necessary. However, some general rules concerning the suppression 
of decoherence can be extracted directly from Q. To avoid decoherence for a 
set of initial states {|m >} we must put the transition probability from |m > 
to all other states equal to zero. This is the case if the energy conservation 
En = cannot be fulfilled or the transitions from any initial |m > to all 
other states are forbidden, i.e. < n|V^|m >= 0. 

It will be useful to notice that the transition probability Pnm can be 
written as an integral 

1 

Pnm=^J Vnm{t)Vmn{0) dt (4) 

where Km(t) =< n\V{t)\m >= exp{i{En- Em)t/h} < n\V\m > and V{t) = 
^iHot/hy^-iHot/h^ each particular application of we have to average the 
transition probability over certain sets of initial and final quantum states 
and then (j3)) becomes an example of Kubo formula for transport coefficients 
expressed as integrals of autocorrelation functions [TKj . 
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2 Dynamical decoupling methods 



It is known from the classical theory of dynamical systems that a fast periodic 
perturbation can either enforce chaotic behavior in the otherwise regular sys- 
tem or in other cases stabilize the system evolution. The later phenomenon 
has been extended into quantum domain and proposed as a method for re- 
ducing decoherence effects in quantum information processing. A rather 
opposite approach utilizes decrease of decoherence rate with an increasing 
time scale of quantum evolution. This type of behavior is characteristic for 
the systems linearly coupled to quantum fields (photons, phonons, etc.). 



2.1 Open system in weak coupling regime 

Consider an open system S weakly coupled to a large reservoir TZ . The total 
Hamiltonian is given by 

H = Hs + Hn + H,^t, H,,t = S ® R . (5) 

The form S ®R for the interaction Hamiltonian (5 = 5'^- system's operator, 
R = -reservoir's operator) is chosen for simplicity only. Most of the 
results can be easily generalized to the interactions of the form J2 ® Ra- 
We assume a discrete spectrum for Hs and a practically continuous one for 

Hs\k>=ek\k> , Hr\E,^>=E\E,^> . (6) 

with 7 describing additional degeneracy of energy levels. 

We compute P^i - "the transition probability per unit time from the state 
|/ > to the state |A; > of the system S" using the Fermi Golden Rule for 
the perturbation applied to the initial state \l > ^\E,'~f > and to a final state 
\k > ®|£",7' >. Averaging over the initial probability distribution for the 
reservoir states a{E,'y) and over all possible final states of TZ we obtain 



Pki='^\<k\S\l> 
n 
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YljdE'j dEa{E,-f)\ < E',i\R\E,-f > \'^S{ei + E - ek - E') . (7) 



X 
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The formula ((Zj) will be used to discuss two oposite strategies of noise reduc- 
tion: minimal decoherence model and "bang-bang" techniques. 
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2.2 Minimal decoherence model 



For many physical models the linear coupling of iS to a quantum bosonic 
field, or, in other words to a system of quantum harmonic oscillators, is an 
important source of noise. For these models the transition between two states 
|/ > and \k > is accompanied by the emission or absorption of a single boson 
(e.g. photon, phonon, etc.) of energy hu. Therefore, the allowed transitions 
for which | < E,'y\R\E' > p > are only between the states of the 
quantum field which differ by a single boson of frequency u, i.e. \E—E'\ = fvjj. 
The energy conservation leads to the condition 



and hence, the transition probability Pu is proportional to the density of the 
bosonic states n(u;) at oj = le^ — ei\/h. Typically, such a density behaves like 
uj^, r > 0, i.e. tends to zero for low frequencies. This suggests a minimal 
decoherence strategy which demands to use systems with almost degenerated 
energy levels and to apply slow gates producing sufficiently small energy level 
splitting - max \ ek — 

There exist natural restrictions for application of minimal decoherence 
strategy. First of all, applying slow gates means using low frequency exter- 
nal (typically electromagnetic) fields for the system control. Low frequencies 
mean long waves of the corresponding fields, which implies difficulties in in- 
dividual control of localized constituents ("qubits") of the quantum device. 
Secondly, beyond the processes of emission or absorption of energy quanta, 
various scattering processes are always present. Their kinematics is com- 
pletely different. For example, in the case of elastic scattering we can have 
|efc — e;| = 0, but the density of outcoming scattered states, which replaces 
the density of bosons from the previous model, is different from zero. 

2.3 Bang- bang techniques 

These techniques utilize the typical behavior of the matrix elements 

< E,'-f\R\E' ,'-f' >. Practically, for all known models of open systems 

< E, 'y\R\E', 7' >~ if E' — E » Ecut for a certain cut-off energy Ecut- 
Introducing a rapidly varying Hamiltonian Hs{t) we can use a modified 
derivation of (|7j). In the interaction picture with respect to the dynamics 
governed by Hs{t) + Hpi we obtain the following time-dependent effective 
Hamiltonian 



= \ek- ei 




Hintit) = S{t) ® R{t) 



(9) 
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where S{t),R{t) are the solutions of the Heisenberg equations 

j^S{t) = l[Hs{t),S{t)] , jR{t) = '-[Hn,R{t)] . (10) 

We use a modification of the formula Q which takes into account the time- 
dependence of Hs(t). This leads to an additional time- averaging yielding an 
autocorrelation function of the type 

F(t)= hm r f(t + s)f{s)ds. (11) 

Applying this generalization of Q to the Hamiltonian Q and averaging 
over the reservoir states as in ((7j) we obtain the transition probability over 
sufficiently long time interval 

Pki = ^^ j duj j dE' j dEa{E)\ < E' ,-f'\R\E,j > \'^Gkiiuj)S{E-E'+hiu) 



r .7 



= ^Y1 / ^^'/ dEa{E)\ < E',y\R\E,^ > I'GkiiiE' - E)/h) (12) 

7' .7 

where Gki{uj) is the power spectrum of the autocorrelation function 

-| /»T /»oo 

Fklit) = lim — < k\Sit + u)\l >< l\Siu)\k > du= GMiuj)e''^' du . 

(13) 

One should notice that by introducing a time- dependent Hamiltonian Hs{t) 
we modify the energy balance. For a constant Hs only the frequencies 
u = {ek — can appear in the power spectrum Gkiiyj). Now, the power 
spectrum can be altered by adding or extracting the energy quanta hw' sup- 
plied by the varying external field present in Hs{t). 

The main idea of the "bang-bang" technique is to shift the support of the 
power spectrum Gki{oj) beyond the cut-off value Ecut in order to make the 
energy conservation condition E' — E = hw impossible to satisfy with uj 
within this support. In other words the overlap of two functions in the final 
integral ()12j) must be very small [TH] . 

In principle it can be done using a Hamiltonian Hs{t) which contains terms 
rapidly oscillating with the frequency Q » Ecut/^- Take a simple example 
of the Hamiltonian periodic in time Hs{t + t) = Hs{t). For the eigenvectors 
of the Floquet operator 

U(r) = U(T,0), U(T)\k>=exp(-'-^)\k> (14) 

n 
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where ^ 

U{t, s) = Texp(^--^ ^ Hs{u)duj (15) 

we have the following representation of the evolution jT2j 

Uit,0)\k>=exp{--^)\Mt)> ■ (16) 

Here (f)k{t) is a periodic (with a period r) function of t with values in nor- 
malized vectors of the Hilbert space. From the periodicity of (f)k{t) it follows 
that (see ^) 

Fki{t) = exp{--(e/ - ek)t}x 

X-/ < + s)\S\Mt + s) >< Ms)\S\Ms) > ds . (17) 
^ Jo 

and the term given by the integral is periodic in t. As a consequence the 
power spectrum possesses the following structure 

Gkiiuj)= V u^5[uj-Uei-ek)-'^] (18) 

n=— oo 

with weights i^ri > depending on the particular choice of Hs{t). There- 
fore, decreasing the period r below the value h/ Ecut and designing a proper 
shape of Hs(t) we can, in principle, reduce transition probabilities Pm. How- 
ever, one should notice that the similar effect can be achieved, simply, by 
performing very fast gates with the duration time tg << h/Ecut- 

Unfortunately, for most of the physically relevant models the cut-off en- 
ergy provides a maximal energy scale for which the given model of a system- 
reservoir interaction is valid. For example, linear coupling to phonons is 
restricted by the Debye energy Ecut = hcuD, in quantum optics such a cut-off 
is provided by the ionisation energy of atoms while in quantum electrody- 
namics one usually takes E^ut = ^eC^ , where rrie is the electron mass. To 
preserve the consistency of the mathematical model we should not consider 
the frequencies satisfying hj » Ecut- If such frequencies can be physically 
realized then usually different mechanisms of decoherence must be also taken 
into account (e.g. nonlinear, 2-phonon processes in solid state). 

Perhaps, the only example for which the above objections do not apply 
is the leading decoherence mechanism in NMR systems. It is related to the 
fluctuations of the local magnetic field possessing a natural cut-off in time- 
scale and hence energy scale as well. The same holds for the "engineered 
noises" used is some of the experiments on quantum error control [T^ . 
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3 Zeno effect 



In the early discussions of the Zeno effect its physical origin was attributed to 
the frequent von Neumann projective measurements, described by the pro- 
jections Pj Pj = l,PjPk = ^jkPj)', such projections combined with the 
Hamiltonian evolution U{t) = exp{ — {i/h)Ht} produce an effective dynamics 
given by the asymptotic formula 

pit) = J2w,it)pW^it) , (19) 
j 

where 

Wj{t) = lim \PjU{t/n)Pjr = P exp{-{i/h)PjHPd] . (20) 

The main consequence of frequent measurements is the stabilization of a 
quantum state which is trapped into a subspace corresponding to the pro- 
jector Pj. This idea has been used, for example, to describe the suppression 
of decay of an unstable particle when continuously observed. Some authors 
went even further to apply Zeno effect in the case of "observation without 
observer", i.e. for a general open system for which interaction with an envi- 
ronment replaces continuous measurements. However, the predictive power 
of such a theory is very weak. Take a neutron, which is unstable as a free 
particle. Inside a nucleus the interaction with other nucleons can stabilize a 
neutron or not. Similarly, a proton, a stable particle by its own, can become 
unstable when put inside a nucleus. To predict the behavior of a nucleon we 
need to know the detailed energy balance of the total composed system and 
the rough picture of "Zeno effect" for a continuously perturbed ("observed") 
particle is rather misleading. 

In more recent papers the phenomenological picture of wave collapse after 
von Neumann measurement has been replaced by more fundamental models 
of a system strongly interacting with an environment represented either by a 
large reservoir (e.g. macroscopic measuring apparatus) or periodic external 
fields. Our aim is to provide an elementary explanation based on the energy 
conservation principle. 

3.1 Open system in strong coupling regime 

We consider again the model of an open system with the Hamiltonian (jSj) 
but with the assumption that the interaction with the environment is strong. 
The Hilbert space basis of a system {\j >} is now determined by S and not 
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by Hs i.e. 

S = Y,s,P, ,P, = \J><J\. (21) 

j 

We use the notation 

H^^^ = Hr + s,R , Hs = Y, ^jP^ + y (22) 

j 

where < >= 0. We can decompose the total Hamiltonian as 

H = Ho + V, Ho = Y,Pj^i^AR + H^R) (23) 

j 

where the off-diagonal part of Hs, denoted by V, is treated as a small per- 
turbation. 

The eigenvectors and the eigenvalues of H^^ are given by 

H'i'^\E,r,J>=E\E,r,J> . (24) 

We can assume that for any j , E > Eg''^ where Eg'''^ is the ground state 
energy of hI^K As we assume the strong coupling between S and R the 
ground state energies E^-''^ depend strongly on the index "j". We can also 
find the eigenvectors and the eigenvalues of Hq 

Ho\j > 0\E,r,J >= (e, + E)\j > ®\E,r,3 > , E> E^^^ . (25) 

Assuming for simplicity the zero temperature reservoir (or in other words 
initial non-degenerated ground state ^ >), we can compute the transition 
probability Pki averaging over the final states of the environment 

27]- 1"^ 

Pki = ^\<k\V\l>?y, dE\< E, ^\R\Ef > \'6{ek + E - - Ef) . 

ft j 

(26) 

If the following threshold condition is satisfied 

Ef^ + efc > + ei (27) 
then the energy conservation forbids the transition |/ >i— |/c >, i.e. Pki = 
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3.2 A model of Zeno effect 



We discuss now a simplified model describing the experimental test of the 
Zeno effect ^H]- A 3-level atom is subjected to the interaction with a strong 
laser field driving the optical transition between the levels 1, 3 and to a 
radiofrequency field causing Rabi oscillations with a frequency Q between 
the levels 1, 2. The model Hamiltonian can be written as 

H = hwislS >< 3| + >< 2| + |2 >< 1\) + hJ2^k[a^ - Fk] [ak - Fk] 

k 

+gni\l >< 3| + |3 >< 1|) ® Y^^fkak + • (28) 

k 

Here a^, are creation and annihilation operators of the electromagnetic 
field and F^ correspond to a classical field describing the strong pumping 
mechanism of the laser. This classical field drives the quantum electro- 
magnetic field into a new ground state - a coherent state \F > satisfy- 
ing ak\F >= Fk\F >. It is convenient to introduce a new atomic basis 
|± >= -^(|1 > ±|3 >), |2 > and new field operators = ak — F^. 
The Hamiltonian can be now written in a form 



H = Ho + V , i7o = P-®(^i'"^ + e-)+P+®(i^i""^ + e+) + P2®i^^"'"^ (29) 
where P± = |± >< ±|,P2 = |2 >< 2| and 

Hiem) ^ J2 Ukblbk , H^r^ = H^^^^ ± g,, Y^^hbu + fkbi) , (30) 



6± = ^ ± g^s Y^ihFu + fkh) , (31) 

k 

V^=-^(l+><-| + |-><+l) 

+^(1+ >< 2| + I- >< 2| + |2 >< -I + |2 >< +1) . (32) 
The ground state energies of the Hamiltonians H^'^"^^ , H^^^ are given by 

i?, = 0,E(-)=EW = -$:^ (33) 

k k 
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respectively with the corresponding coherent ground states \F >, > 
satisfying 

6fc|F>=0,6,|FW>=TA|F(±)> . (34) 

Assume that gis '^^{fkFk + IkFk) > 0. Then, the energy of the ground state 
|— > > of Hq, for a strong laser field {F^}, is much lower than the 

energy of the Hq- eigenstate |2 > <S)\F >. Therefore, according to ()2(jp27|l 
the transition from the state |— > to the state |2 > driven by the Rabi term 

>< 2| + |2 >< 1|) is forbidden by the energy conservation principle. 
This phenomenon is interpreted as the Zeno effect caused by continuous 
measurements performed on the levels 1 and 3 by means of the interaction 
with a strong electromagnetic (laser) field. 

The main difficulty in application of the Zeno effect in quantum infor- 
mation processing is the presence of a strong coupling with environment. 
The structure of the unperturbed Hamiltonian ()23I25|) suggests that a sub- 
system S becomes a dressed system and the very definition of the collection 
of "individual qubits" with its own controlled unitary dynamics becomes 
questionable. 

4 Decoherence free subspaces and subsystems 

Consider a model of an open system S weakly interacting with a reservoir R as 
in Section 2.1. The physical mechanism leading to decoherence-free subspaces 
(DFS) is entirely due to the vanishing matrix elements < k\S\l >= for all 
pairs of states, such that \l > belongs to a certain subspace TC^^^ of the 
Hilbert space of S and |/c > is an arbitrary state orthogonal to |/ >. This 
happens if and only if S\l >= s\l > for all \l >G Ti.^^^ , i.e. the operator S 
acts as a scalar on the subspace 7i^^^ . In such a case we can say that the 
transitions |/ \k > are forbidden. 

4.1 Symmetries and forbidden transitions 

Typically, decoherence free subspaces are associated with the symmetries 
characterizing the system and its interaction with the environment jSlEJlE]- 
Assume that these symmetries are given in terms of a generally reducible 
representation of a certain semi-simple Lie group acting on the Hilbert space 
of the system S. This representation is given by a basis of operators = 
generating the representation of the corresponding Lie algebra and choosen 
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in such a way that the Casimir operator C can be written as [20] 

C = E^J- (35) 

Denote by Hq the subspace spanned by the eigenvectors of C with the eigen- 
value zero. It follows that for any (p G Tio we have 

=< 0|C|0 >= J2 ^X^<p = 0. (36) 

Therefore for any interaction Hamiltonian (jSJ which is of the form 

H,^t = J2x^®R^ (37) 

the subspace T^o is decoherence-free. 

Typical examples are provided by systems with rotational (or isospin) sym- 
metry. They are equipped with a (generally reducible) representation of the 
SU{2) group and corresponding "angular momentum" operators [Jk,Ji] = 
iekimJm, k,l,m = 1,2,3, and the Casimir operator = + J| -|- J| with 
the eigenvalues j(j + 1), j = 0, 1/2, 1, .... In this case all singlet states corre- 
sponding to j = span the DFS for the interactions with an environment of 
the form Y,k Jk® Rk- 

4.2 Collective decoherence 

The most popular model leading to DFS consists of 2N qubits coupled to en- 
vironment by means of a collective interaction ( for the related superradiance 
and subradiance phenomena, see [2III2Z|) 

2N 3 

Hint = Y.ll''k^®^^ (38) 
i=l k=l 

where a^^\ are Pauli matrices for the j-th qubit. The operators Jk = 
\ Sj=i '^k \ ^ = 1, 2, 3 generate a reducible representation of SU{2) an there- 
fore the 2A^-qubit singlet states span DFS of dimension {2N\)/{N + l)!iV!. 
Here the mechanism leading to DSF is due to the symmetry of the interaction 
Hamiltonian with respect to qubits permutations. 
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4.3 Decoherence-free subsystems 



A generalization of DFS leads to the notion of decoherence-free subsystemsf^ . 
One assumes that the operator S can be represented in the following block- 
diagonal form 

S = ^Inj^Sj. (39) 

J 

corresponding to the decomposition of the Hilbert space 

= C"'^ ® C^' (40) 
J 

typically associated with reducible representations of compact Lie groups. It 
means that 5* acts as a scalar on Hilbert spaces C"-^ which now correspond 
to certain fictitious subsystems. According to the previous discussion these 
subsystems are decoherence-free. 

As argued in the idea of topologically protected states [2^ is also 
closely related to decoherence-free subsystems. 

4.4 Physical implementation 

The engineering of the interaction Hamiltonian in real systems is a very 
difficult task. Namely, the permutation invariance can be only approximative 
and the operators Rk always depend on j also. To reduce this dependence one 
should place all qubits within a space region of the linear dimension which is 
small in comparison with the typical wavelength of the reservoirs fluctuations 
coupled to the qubit system. On the other hand, a dense packing of qubits 
leads to unwanted interactions between them which are not permutation 
invariant and hence spoil the effect of collective decoherence suppression. 

5 Conclusions 

We have shown using the elementary quantum mechanical time-dependent 
perturbation theory that the different methods of controlling decoherence 
are based on the combination of the energy conservation principle with the 
concept of forbidden transitions usually related to symmetries of the system. 
All those methods possess their own rather restricted ranges of applicability 
and demand often contradictory means to be implemented. For instance, the 
Zeno effect is based on a strong interaction with an environment while the 
other techniques rely on the weak coupling assumption. Minimal decoherence 
technique demands coupling to low frequencies of the reservoir, similarly to 
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DFS based on permutation invariance guaranteed by the coupling to long 
wave fluctuations. On the other hand "bang-bang" techniques employ high 
frequency regime. One can expect that only an optimized combination of 
different techniques might be, to some extent, successful in noise reduction 
for microscopic quantum devices (see [22] for the overview of hybrid methods). 

The author is grateful to Daniel Lidar and Marco Piani for very useful 
comments. This work was supported by the Polish Ministry of Scientific 
Research under Grant 2P03B 084 25 and by EC grant RESQ IST-2001-37559. 
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